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PACS 73 . 20 . Fz - Weak or Anderson localisation 

PACS 02.50.-r - Probability theory, stochastic processes, and statistics 

Abstract. - A one-dimensional model of classical diffusion in a random force field with a weak 
concentration p of absorbers is studied. The force field is taken as a Gaussian white noise with 
(4>{x)) = and {4>{x)4>{x')) — gS(x — x'). Our analysis relies on the relation between the Fokker- 
Planck operator and a quantum Hamiltonian in which absorption leads to breaking of supersym- 
metry. Using a Lifshits argument, it is shown that the average return probability is a power law 
(P(x,t\x,0)} ~ t~\/ 2 P/9 (to be compared with the usual Lifshits exponential decay exp — {p 2 t) 1 ^' i 
in the absence of the random force field). The localisation properties of the underlying quantum 
Hamiltonian are discussed as well. 
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Introduction. — Classical diffusion in a random force 
field is encountered in several physical contexts of statis- 
tical physics related to anomalous diffusion or glassy dy- 
namics due to the presence of disorder. As introduced by 
Sinai in the seminal work [1], it may be modelised via a 
Langevin equation x(t) = 2(j>(x(i)) + V2rj(t), where <f>(x) 
is a quenched random force field with short-range corre- 
lations, and rj(t) a Langevin force (a normalised Gaus- 
sian white noise). By now its large-time properties in 
one dimension are well understood : in the absence of 
a global drift, the random force leads to anomalous diffu- 
sion characterised by the scaling of the distance with time 
x(t) ~ In 2 t. Normal diffusion properties are only recov- 
ered for a sufficiently large drift, however an intermediate 
regime reveals several interesting phases (see Rcf. [2] for 
a review). Many approaches to this problem have been 
developed : a probabilistic method [2] (a continuous ver- 
sion of the Dyson-Schmidt method [3]), Berezinskii dia- 
grammatic techniques [4, 5] as well the replica method [2] . 
Moreover more recently, interesting features of this model 
like aging properties were analyzed by means of Ma- 
Dasgupta real-space renormalisation group methods [6] . 



In this letter, we extend the analysis to random me- 
dia containing randomly spread absorbers. Our approach 
relies on the Fokker-Planck equation (FPE) 
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describing classical diffusion in a force field (j>(x) in the 
presence of an absorber density A(x). P denotes the (con- 
ditional) probability P(x, t\x' , 0) to find a particle at x at 
time t which has started from x' at t! = 0. We consider 
the random force field to be 4>{x) a Gaussian white noise 
of zero average (4>(x)) = and (4>{x)4>(x' )) = g5(x — x') 
(throughout the paper, (• • •) will denote averaging with re- 
spect to the quenched disorder : force field and absorbers). 
A(x) describes absorbers at locations x n with annihilation 
rates a n > 0, independently and uniformly distributed for 
a concentration p ; thus we write A(x) = J2n a " S(x — x n ). 
The effects of the random force field and the random local 
annihilation rates are well known when described sepa- 
rately. Let us first review known results for the average 
return probability (x = x'). 

(i) g — & a. — : In the absence of random force 
fields and impurities, the conditional probability reads 
P(x,t\x,0) = -J=. 
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(ii) g ^ & a. — : For classical diffusion in a random 
force field (Sinai problem) [1,2,6], the decay is 

(P(x,t\x,0)) ~ . 9 ln-V*), (2) 

r — >co 

and hence much slower than This behaviour is 

related to the aforementioned typical distance x(t) ~ 

(iii) g = a ^ : For free diffusion with a weak 
absorber concentration we have 

(P(x,t\x,0)) ~ p(p 2 1/6 e- 3 ^) 2/3 (" 2t ) 1/3 . (3) 

i — >oo 

The rapid (exponential) decay is mostly explained by the 
decay of the survival probability J dx (P(x, t\x' , 0)) ~ 
exp — (p 2 t) 1 / 3 (probability is not conserved in the pres- 
ence of absorption) (see [7] for a review on Lifshits tails). 
As switched on from free diffusion, the random force field 
strongly slows down the decay of the probability from 
\j\fi to 1/ln 2 t, whereas absorbers tend to accelerate the 
decay from a power law to exponential exp — t 1 / 3 . The aim 
of the present paper is to study the interplay between the 
random force field and randomly dropped absorbers. 

From FPE to Schrodinger equation. Our anal- 
ysis relies on the well-known relation between the FPE (1) 
and the Schrodinger equation —dttp = Hip for : 

H = --^2 + 0(x) 2 + <j>'{x) + A(x) = H snsy + A(x) . (4) 
dor 

A mapping between the two equations is constructed 
via a non-unitary isospectral transformation P(x, t) = 
ip(x,t)exp J (p(x')dx'. In the case A = 0, this leads 
to the Hamiltonian ff SU8y = + <!>)(-£ + <j>)- This 
factorised, supersymmetric, structure is responsible for a 
positive spectrum. When ^ is a white noise of zero mean, 
the spectrum presents a Dyson singularity at the band 
edge (E = 0) [2,4,5,8]. The presence of the absorption 
A(x) in Hamiltonian (4) breaks the supersymmetry. Spec- 
tral and localisation properties of H were investigated in 
Rcf. [9] for the case of A(x) and <p(x) both being white 
noises. This work focused on the mechanism leading to 
a (stochastic) supersymmetry breaking, giving rise to a 
lifting of the Dyson singularity and a breaking of the de- 
localisation transition at energy E = 0. This first study is 
connected to the high density limit p > a„ of the model 
studied in the present paper. However Ref. [9] has led to 
the conclusion that the low density limit p <C a n studied 
below is more relevant in the context of classical diffusion. 

In order to study diffusion properties via (4) , let us recall 
that we may relate the return probability, averaged over 
the realisations of the random functions 4>(x) and A(x), 
to the Laplace transform of its density of states (DoS) : 

(P(x,t\x,0)) = J dE p(E) e~ Et . (5) 

Having this relation in mind, we now construct a Lifshits 
argument for the DoS p(E) of the Hamiltonian (4). 



Free diffusion with random absorbers (g = & a ^ 0). 

We first recall the famous Lifshits argument [3, 10] in 
the absence of the random force field = 0. Low-energy 
states are due to the formation of large impurity-free re- 
gions. Let us denote by £ the distance separating two 
neighboring impurities. For l/£ ~ p <C a n they im- 
pose on the wave function to vanish at their location 
(indeed this holds rigorously for a n — » oo) and the in- 
terval gives rise to a low-energy state E\ ~ (ir/l) 2 . In 
terms of classical diffusion, the survival probability of a 
diffusive particle released in such an interval decays as 
e -7r-t/f _ Hence the probability for a low-energy state is 
related to the probability of a formation of a large interval : 
Proba[£i < E] = Proba[£ > tt/VE}. Since the distribu- 
tion of £ is pe~ pl we recover the Lifshits singularity [3,10] 
for the integrated density of states (IDoS) per unit length 
N(E) ~ q-kp/^/e _ Using a steepest descent method, we 
can relate this low energy behaviour to the large time be- 
haviour of (P(x, t\x, 0)) and recover eq. (3) 1 . 

Let us work out a more precise argument : in the limit 
of high impurity weights, a n — * oo, an interval of length £i 
between two impurities yields a contribution Afo{E;£i) = 
XT=i e i E - to the IDoS : 

where the average is taken with respect to I (this approx- 
imation corresponds to the "pieces model" of Ref. [3,12]). 
We have N{E) = p £~ x ^/Ve de P e ~ P ' what y iclds I 13 ] 

N(E) ~ -4= for s/E, p < a. (7) 

Random force field with absorbers (j / S a / Oj. 

We now apply the same argument to the Hamiltonian 
(4) in order to obtain the low energy DoS. Due to the 
supersymmetric potential (j){x) 2 + <fi'(x), the energy levels 
E n associated to an interval of length £ differ from (mr/£) 2 , 
and rather are distributed according to some nontrivial 
laws W n (E;£) obtained in Ref. [14]. Similarly to (6), we 
have 

N(E) ~ p(Af S usy(E;£)) e , (8) 

where Af susy (E; £) is the IDoS of H susy on [0,£] 
for Dirichlet boundary conditions. We use the 
decomposition over the distributions of eigenvalues 
Kn Sy (E;£) = Y.n^Jv AE ' W n{E';£). Following [14], 
in the limit g£ 3> 1 and for E -C <? 2 , we may 
write the distribution of the n-th energy level as 
W n (E;£) ~£N^ usy (E) w n (i N susy (E)), where iV susy (£) = 

x This picture can be generalised in higher dimensions where the 
main exponential behaviour N(E) ~ Q~ pE is due to low lying 
states of energy E ~ 1/L 2 in large regions of volume L d free of 
impurity associated to probability Q~P L [3,10]). This leads to 

2 d 

(P(x,t\x,0)) ~ exp — p d + 2 t d + 2 . The preexponential factor of the 
IDoS has been studied by instanton techniques [11]. 
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limp 



W S usy(E;£) 



[2, 8] is the 
; Jq(x) and 



SUSy ^' V ~ Jo(-^E/9) 2 +N (VE/ g )i 

IDoS per unit length for an infinite volume 
Nq(x) are the Bessel functions of first and second kind, 
respectively. Contrary to Af SUS y {E; £) , the IDoS per unit 
length N susy (E) is insensitive to boundary effects. Finally 



N(E) 



f'N. 



uay(-E) 

dx w n {x) 



(9) 



We use the integral representation [14] vu n (x) = 
f n — tin. I— , where B is a Bromwich contour (axis go- 
ing from c — ioo to c + ioo with all singularities of the 
integrand at its left), and obtain : 



0.30 



0.10 - 




N(E) 







n=l 



2i7T qcosh^^/q 



(10) 



If we permute order of integrations, pcrfom the integral 
with respect to £ and the summation we find 



N(E) ~ p 



1 



1 



dq 

B 2i7r p/N susy (E) - q sinh 2 y/q 



(11) 



Notice that this procedure only converges for a Bromwich 
contour with < c < q* = p/N susy (E), what always can 
be attained via a suitable contour deformation. Applying 
the residue theorem to the simple pole at we obtain 
our main result for the low energy IDoS per unit length 



N(E) 



sinh 2 y/p/N BVSy (E) 



for E <C g , a 
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The aforementioned condition g£ 3> 1, where £ denotes 
the length of an interval, turns out to be a low density 
condition p <C g under which (12) is valid. For complete- 
ness, note that E 3> g 2 , a 2 corresponds to the pertur- 
bativc regime where we recover the free IDoS N(E) ~ 
—y/~E. Eq. (12) allows to identify the energy scale E c — 

-\/2g/p SC p ara ting two regimes. In the intermediate 



g*e~ 

encrgy range, E c <C E <C g 2 , we recover the IDoS of H s 



N(E) ~ N susy (E) 



2g 



ln 2 (g 2 /E) 



(13) 



It is only in the narrow region < E < E c that the scalar 
potential A(x) affects the spectrum, for E <C E c : 



N(E) ~ 4p 



(14) 



Interestingly, we notice that this power law behaviour 
stems from the ground state energy E\ of each interval 
between consecutive impurities. Indeed, we may check 
that eq. (14) can be obtained from N(E) ~ p(8(E — 
Ei[<p(x), £]))</,. £, by means of the distribution Wi(E;£), in- 
volving Wl (a;) = YZ=^T + ^^ ' x M- & 



Figure 1: Comparison between numerical results (triangles, di- 
amonds & stars) and analytical results (lines). Inset : IDoS 
for g = 1 with p ranging from 0.01 to 0.5. 

priori it is far from being obvious that the analysis can 
be restricted to the ground state of each interval since the 
distributions W n (E; £) have strong overlaps [14] . 

Again, combining (14) and a steepest descent argument 
for (5), we relate the power law behaviour of the IDoS to a 



power law decay of the particle density for t ^> p-e"^ 



(P(x,t -> oo|x,0)) 



4 P r( v /277^+l) 
(gH)VW 3 



(15) 



The crossover time scale t c = g~ 2 exp y / 2g/p corresponds 
to the time needed by the random particle released in the 
random force field to reach the nearest absorber : x(t c ) ~ 
\jp where x(t) ~ g~ 1 \n 2 (g 2 t). The return probability 
(15) decays faster than (2) but slower than the exponential 
decay (3) in the absence of the random force field. This 
provides the answer to our initial question. 

Localisation. — We now analyse the localisation 
properties of the underlying quantum Hamiltonian (4). 

(i) g & a. — : In the absence of absorption, A(x) = 
0, the Lyapunov exponent (inverse localisation length) 
is exactly known [2] 7 SUS y(-E) = — f fcgx In [Jo (k/g) + 
reaches a finite value at high energy 
g/2 and vanishes at zero energy : 

9 



N 2 (k/g)}. It 

7susy(£> g 2 ) = 



7susy (S^0) 



Hg/k) 







(16) 



(m) g = & a ^ : On the other hand, in the absence 
of supersymmetric noise, <j)(x) = 0, and for a low density 
of impurities, the Lyapunov exponent is given by 

7scalar(£) ^ § ln[l + (f ) 2 ] for k » p [3,15]. At high 
energy it leads to the well known linear increase of the 
localisation length as a function of the energy [16] l/£\ oc = 
7scaiar(£') ~ pa 2 / (8E) for E — > oo. At zero energy it 
reaches a finite value given by [3, 17] 

(E = 0) ~ p [\n{a/p) - C] for p < a , (17) 
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where C ~ 0.577 is the Euler constant. 

(Hi) g y£ & a ^ : We now turn to the case where <p 
and A both differ from zero. The high energy (perturba- 

2 

tive) expression is : j(E) ~ + 1 . It corresponds to the 
addition of the pcrturbative expressions for 7s Ca iar(£') and 
7susy(^) (a similar result was obtained in Ref. [9] when A 
is a Gaussian white noise). We can see on figure 2 that, for 
a fixed g, 'y(E) (obtained numerically) slowly converges to 

7susy 

(E) as p is decreased. The introduction of the scalar 
potential A(x) breaks the dclocalisation at E = obtained 
for the supersymmctric Hamiltonian. 



<« <<<<f=3<<i <<<<<] <<<] < ] < . 
p=0.2 



>£0.6 




Figure 2: Low energy Lyapunov exponent as a function of en- 
ergy for p ranging from 0.2 to 0.01 ; with g = 1. Black line : 
7suay(-E). Green dashed line : perturbative result for p = 0.2. 

Let us now turn to the detailed analysis of the zero 
energy Lyapunov exponent. For this purpose it is con- 
venient to convert the Schrodingcr equation Hip = Eip 
into a stochastic differential equation for the Ricatti vari- 



ip' /ip — (p : this latter obeys 4-z 



-E 



able z 

2z<p[x) + A(x). This Langevin like equation can be re- 
lated to a Fokker-Planck type equation for the distribution 
d x T{z-x) = d z [(E + z 2 )T(z;x)} + 2gd z [zd z [zT(z; x)}] + 
p[T(z — a; x) — T(z;x)], where the first term is a drift 
term related to the force field F(z) = —E — z 2 , the second 
term a diffusive term and the last term is a jump term 
originating from the scalar impurities A(x). For E > 
the Ricatti variable is driven from +oo to — oo in a finite 
"time" x. The steady current of z(x) corresponds to its 
number of divergences per unit length, i.e. the number of 
nodes of the wave function ip per unit length, that is to 
N(E). The distribution reaches a stationary distribution 
T(z; x — ► oo) = T(z) for a steady current —N(E) : 

N(E) = (E + z 2 )T(z) + 2gz^-[zT(z)} — p f dz'T(z') 

dz J z - a 

(18) 

The IDoS is given by normalising the solution of this 
integral equation. Given T(z), the IDoS can be ex- 
tracted from the distribution thanks to the Rice formula 
z 2 T(z). Lyapunov exponent is given 



N(E) = lim*-, 

by [3] y(E) = (z) = lim^co J^ dzzT(z) 



For E = and positive jumps a > 0, the Ricatti is 
constrained to belong to R + since both the "force" F(z) 
and the effect of the multiplicative noise <fi(x) vanish for 
z = : we recover N(E = 0) = 0. In the low density limit 
p — > the Ricatti variable is driven to z ~ and eventually 
reinjected at z ~ a with a "rate" p. We can write that 
the negative current due to the "force" F(z) and the noise 
4>{x) is equilibrated by the positive current p due to the 
jumps : p ~ z 2 T(z) + 2gz-^[zT(z)] for z G [z c ,a]. We 
have introduced a cutoff z c where the jump occur almost 
surely (the force F(z) and the multiplicative noise vanish 
effectively as z — > 0). The solution of this equation is 



T(z) ~ o -T g - Z ^ zl29 f Zr ip-e z /2g - This distribution may 
be approximated by T(z) ~ -^\ii(z/z c ) for z c < z < g 
and T(z) ~ for g < z < a (and T(z) ~ elsewhere). 



Normalisation gives In (g / z c 



-f. Using 7 



7(0) 



pg + pln(a/g) 



we get 



(19) 



The second term is reminiscent of the result obtained in 
the absence of the supersymmctric noise. However the 
dominant term involves a nontrivial combination of p and 
g. In the region of rarefaction of eigenstates (E ~ 0), 
the localisation length reads i\ oc = I/7 ~ j y/p/s 1/P- 
The scalar impurities breaks the delocalisation transition 
of the supersymmctric Hamiltonian, but quite surprisingly 
the localisation length is much smaller than the inverse 
density of impurities. 

Numerics. — We now turn to a numerical analy- 
sis in order to check our analytical results for spectrum 
and localisation and explore more precisely their valid- 
ity range. We analyse the Hamiltonian (4) for 4>(x) = 
J2n^nS(x-x n ) and A(x) =Yl n ot n S(x — x n ). This choice 
of modelisation of the random force field <fi allows to deal 
with a continuous description. The locations x n are uni- 
formly distributed and uncorrelated, with densities p§ for 
those of <p(x) and p for impurities of A(x). The precise 
shape of the distribution for the dimensionless weights 
A„ is not important, as long as it satisfies (A„) = 
and (A„) finite. We choose a symmetric exponential law 
p(X n ) = 7^e~' A "'/ A . The process cf>(x) behaves as a Gaus- 
sian white noise in the limit — > 00 and A — > 0, with 
9 = P4>( X l) = 2 P^ 2 nxe d. 

Phase formalism. Let us now explain how we can ob- 
tain the spectral density from the phase formalism [3,16]. 
In the equation Hip — Eip, we replace the couple of vari- 
ables (ip,ip') by the couple (9, £) defined as ip = e^sin^ 
and ip' — (pip = ke^ cos 9, where E = k 2 . The two variables 
obey the differential equations [9] : 



d0 
d.r 
d£ 

d.r 



A(x) 



2k 



sin 26> 



+ <p(x) sin 29 
>(x) cos 29 



(20) 
(21) 



We solve these equations as follows 
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Figure 3: Dashed lines correspond to N(E) oc E^/ 2p ^ a . Impu- 
rity weight is a — 50 and the number of impurities is Ni = 10 7 
(Ni = 10 s for the largest densities and smallest energies). The 
red line is eg. (7) and the red stars are the corresponding nu- 
merical calculation, for p = 10. 

• We denote t n = x n +i — x„. The evolution of the 
variables on an interval free of impurity is — 
9+ = k£ n and — £+ = 0, where 0± = 0(x„) and 
£^ = are values just before (— ) or right after 
(+) the n-th impurity. The size of these intervals 
is distributed according to the Poisson law P(t) = 
Ptote~ Ptot£ for a density p tot = p + p . 

• The effect of a (5-peak of A(x) is given by cotg6>+ — 
cotgfl" = ^ and £+ - C = ln|^ = ± In [l + 



^ sin 29- + % 

The effect of 
tan#+ = 
iln [e 2A 



sin 2 9~] (see Ref. [17] 



tan 9~ e 2A 



a (5-peak of 
and £+ - 



is given 

1 l sin20~ 

2 111 si „2fl + 



by 



sin 



cos 2 *?" (see Ref. [15]). 



It is worth noticing that the two effects on the envelope of 
the wave function are roughly given by A£ n = — £~ ~ 
lna„ and A£ n — |A„|. 

The IDoS corresponds to the number of zeros of the 
wavefunction of energy E, i.e. the number of times the 
cumulative phase 0{x) coincides with an integer multiple 
of 7r. Therefore, a convenient way to get the IDoS nu- 
merically is [3, 16] N(E) = lim-c^oo Moreover, the 
damping of the envelope is characterised by the Lyapunov 
exponent 'y(E) = lim x _ s . 00 providing a definition of 
the inverse localisation length (£\ oc = 1/7) . 

Ricatti variable. The equations of the phase formal- 
ism are singular in the limit E — > 0. In order to ob- 
tain the zero energy Lyapunov exponent, it is more sim- 
ple to perform the analysis in term of the Ricatti vari- 
able. Let us denote z~ = z(x^) its values before/after 
the impurity n. Between two impurities the evolution is 
given by arctan (z~ +1 /k) — arctan (z+/fc) = —k£ n , i.e. 
l/^n+i — V z n = for E = 0. Through a scalar impurity 



of A(x) we have z„ — z~ = a„ and through an impurity 
of (j){x) we have z+ = z~ e~ 2A " . IDoS and Lyapunov ex- 
ponent can be extracted from the stationary distribution 
of the Ricatti as explained above. 

Results. As a first check, we consider the case with 
= 0: we compare numerics for A = to the Lif- 
shits singularity (7) (red stars and red line on figure 1) 
and obtain a good agreement (some deviations appear 

for larger k since N(E 3> a 2 ) ~ ^). Next, we treat 
the case with A = and check the numerics (black tri- 
angles on figure 1) against the analytical expression re- 
called above [2,8] (black continuous line) N susy (E — > 0) = 

The agreement to the theory 

is 



3/2 | r>( 



k 2 ^ 

In' 2 k i 



[ln(2 9 /fc) 

shows that our modclisation of Gaussian white noise 



adequate. 



I. 



I 



pig 



Figure 4: Numerical result for function B(x). Lines are 
indicative (they correspond to er Ax and e~ 2x /A). 



Finally, we combine both noises (green diamonds on 
Fig. 1 and inset), and check against the power law N(E — ► 

0) oc E^/ 2p ^ 9 : the exponent fits very well within a surpris- 
ingly wide range (see figure 3 for densities ranging from 
p = 0.01 to 5). We insist that slopes in the log-log plot 
are not fitted but directly compared to \j2pj g (straight 
dashed lines). We however observe that, apart for the 
lowest densities, the prefactor significantly differs from 4p, 
eq. (14). Since the IDoS reaches a finite limit for a — » oo, 
the additional dimensionless factor is a function of the ra- 
tio p/g only. We conclude that numerics suggests the form 

N(E -+ 0) ~ 4 P B(p/g) («)\^ with B(x - 0) = 1. 
The function B(x) is extracted from numerics and plotted 
on figure 4. For p/g = 10, the numerical precision does 
not allow a convincing fit for the lowest energies (green 
diamonds of Fig. 3). We have however checked that the 
noise (f>(x) still affects the IDoS which has not yet reached 
the Lifshits result (7) (red stars and line). 

The energy dependence of the Lyapunov exponent is ob- 
tained from the phase formalism. The results are plotted 
for different densities p on figure 2. In a second step we 
analyze the dynamic of the Ricatti variable for E = 0. The 
stationary distribution is plotted on the inset of figure 5 on 
which we can check the crossover between the behaviours 
T{z) oc 1/z and T(z) oc 1/z 2 occuring for z ~ g = 1 
(note however that prefactors do not fit with the one de- 
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rived above). The distribution is used to compute the zero 
energy Lyapunov exponent. The analytical result (19) is 
compared to the numerical result and satisfactory coin- 
cides (Fig. 5). 




icr 2 icr 1 10° -io 1 

P 



Figure 5: Lyapunov exponent at E — 0. Green dashed line : 
in the absence of super symmetric noise, for a = 50. Numerics 
(circles) is compared to equation (17). Red : with cj){x) (g = 1). 
Comparison with expression (19). Inset : Stationary distribu- 
tion of the Ricatti variable in this latter case (p = 0.05) ; dashed 
blue line is oc 1/z and red line oc 1/z 2 . 

Conclusion. — We have studied the average return 
probability for one-dimensional classical diffusion in a ran- 
dom force field and the presence of absorbers at weak con- 
centration (p — ► 0), yet with strong absorption rates. We 
have shown that absorption only takes place above a very 

large time scale t c = g~ 2 ^J 2g l 9 . The well-known Sinai 
decay (P(x,t\x,0)) ~ 2g\n~ 2 (g 2 t) holds for t < t c and is 

replaced by the power law (P(x,t\x,0)) ~ 4p (g 2 t)~^/^> 
for t > t c . Whereas a simple guess would have been to put 
the crossover between (3) and (15) at g ~ p, the power law 
(15) persists numerically up to large ratio p/g. It would 
be interesting to understand more carefully the crossover 
and obtain analytically the prefactor B(p/g) not predicted 
in our calculation. Another interesting issue would be 
to investigate the fluctuations of the return probability 
P(x,t\x,0) over disorder configurations ; such a question 
is related to the characterisation of the fluctuations of the 
local DoS of the quantum Hamiltonian, a question studied 
for high energies for the scalar noise in Ref. [18] and for 
the supersymmctric noise in Ref. [19]. 

The localisation properties of the underlying quantum 
Hamiltonian have been considered, too. In particular, the 
scalar impurities lifts the divergence of the localisation 
length (inverse Lyapunov exponent) at energy E = and 
leads to a localisation length £\ oc ~ 1/^fgp for p — > 0. 
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